On semidefinite programming relaxations of the 
traveling salesman problem 

Etienne de Klerk* Dmitrii V. Pasechnik"'" Renata Sotirov* 

February 11, 2009 



We consider a new semidefinite programming (SDP) relaxation of the 
symmetric traveling salesman problem (TSP), that may be obtained via 
an SDP relaxation of the more general quadratic assignment problem 
(QAP). We show that the new relaxation dominates the one in the paper: 

[D. Cvetkovic, M. Cangalovic and V. Kovacevic-Vujcic. Semidefi- 
nite Programming Methods for the Symmetric Traveling Salesman Prob- 
lem. In Proceedings of the 7th International IPCO Conference on Integer 
Programming and Combinatorial Optimization, 1999, 126-136, Springer- 
Verlag, London, UK.] 

Unlike the bound of Cvetkovic et al., the new SDP bound is not dom- 
inated by the Held-Karp linear programming bound, or vice versa. 
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1 Introduction 

The quadratic assignment problem (QAP) may be stated in the following form: 



where A and B are given symmetric n x n matrices, and LT n is the set ofnxn 
permutation matrices. 
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It is well-known that the QAP contains the symmetric traveling salesman 
problem (TSP) as a special case. To show this, we denote the complete graph 
on n vertices with edge lengths (weights) Dij — Dji > (i ^ j), by K n (D), 
where D is called the matrix of edge lengths (weights). The TSP is to find 
a Hamiltonian circuit of minimum length in K n {D). The n vertices are often 
called cities, and the Hamiltonian circuit of minimum length the optimal tour. 

To see that TSP is a special case of QAP, let C% denote the adjacency matrix 
of C„ (the standard circuit on n vertices) : 



Ci 



1 

1 

1 





1 







1 
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Now the TSP problem is obtained from the QAP problem fl} by setting A 
and B = C\. To see this, note that every Hamiltonian circuit in a complete 
graph has adjacency matrix XC\X T for some X E Ii n . Thus we may concisely 
state the TSP as 



TSP, 



opt 



mm trace 
xen n 



X -BXC X X T 



(2) 



The symmetric TSP is NP-hard in the strong sense [2D], and therefore so is 
the more general QAP. In the special case where the distance function of the 
TSP instance satisfies the triangle inequality (metric TSP), there is a celebrated 
3/2-approximation algorithm due to Christofides [DJ. It is a long-standing (since 
1975) open problem to improve on the 3/2 constant, since the strongest negative 
result is that a (l+l/219)-approximation algorithm is not possible, unless P=NP 

In the case when the distances are Euclidean in fixed dimension (the so-called 
planar or geometric TSP), the problem allows a polynomial-time approximation 
scheme pQ. A recent survey of the TSP is given by Schrijver [24], Chapter 58. 



Main results and outline of this paper 

In this paper we will consider semidcfinite programming (SDP) relaxations of 
the TSP. We will introduce a new SDP relaxation of TSP in Section [2] that 
is motivated by the theory of association schemes. Subsequently, we will show 
in Section [3] that the new SDP relaxation coincides with the SDP relaxation 
for QAP introduced in [25] when applied to the QAP reformulation of TSP in 
(J2J . Then we will show in Section [J] that the new SDP relaxation dominates 
the relaxation due to Cvetkovic et al. [5]. The relaxation of Cvetkovic et al. is 
known to be dominated by the Held-Karp linear programming bound [6] [15] , 
but we show in Section [5] that the new SDP bound is not dominated by the 
Held-Karp bound (or vice versa). 
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Notation 

The space of p x q real matrices is denoted by W xq , the space of k x k symmetric 
matrices is denoted by , and the space of kxk symmetric positive semidefinite 
matrices by S£ ■ We will sometimes also use the notation X y instead of 
X e <Sjjr, if the order of the matrix is clear from the context. By diag(A) we 
mean the n- vector composed of the diagonal entries of X € <S n . 

We use /„ to denote the identity matrix of order n. Similarly, J„ and e n 
denote the n x n all-ones matrix and all ones n- vector respectively, and 0„ X n is 
the zero matrix of order n. We will omit the subscript if the order is clear from 
the context. 

The Kronecker product A® B of matrices A £ W xq and B S ]R rxs is defined 
as the pr x qs matrix composed of pq blocks of size r x s, with block ij given 
by AijB (i = l,... ,p), (j = l,...,q). 

The Hadamard (component-wise) product of matrices A and B of the same 
size will be denoted by A o B. 



2 A new SDP relaxation of TSP 



In this section we show that the optimal value of the following semidefinite 
program provides a lower bound on the length TSP op t of an optimal tour: 



mm 
subject to 



| trace (DX^) 
> 0, 

EUx (k) = J-i, 

I + J2 d k =iCos(^)X^h0, 



i,...,d 



2=1, 

k = 1, 



,d 
, d, 



(3) 



where d = \ hn\ is the diameter of C n . 

Note that this problem involves nonnegative matrix variables X^\ . . . , X^ 
of order n. The matrix variables X^ have an interesting interpretation in terms 
of association schemes. 



Association schemes 

We will give a brief overview of this topic; for an introduction to association 
schemes, see Chapter 12 in [TP] , and in the context of SDP, [TTj . 

Definition 1 (Asssociation scheme). Assume that a given set ofnxn matrices 
Bo, . . . , B t has the following properties: 

(1) Bi is a — 1 matrix for all i and Bq = I; 

(2) EiBi = J; 

(3) Bi = B% for some i*; 
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(4) BiBj = B 3 Bi for all 

(5) B^ espan{Bi,...,B t }. 

Then we refer to {-Bi, . . . , B t } as an association scheme. If the B^s are also 
symmetric, then we speak of a symmetric association scheme. 

Note that item (4) (commutativity) implies that the matrices B\ , . . . , B t 
share a common set of eigenvectors, and therefore can be simultaneously diag- 
onalized. Note also that an association scheme is a basis of a matrix-* algebra 
(viewed as a vector space). Moreover, one clearly has 

trace(SiBj) = if i ^ j. 

Since the B^s share a system of eigenvectors, there is a natural ordering of their 
eigenvalues with respect to any fixed ordering of the eigenvectors. Thus the last 
equality may be interpreted as: 

^A fe (B i )A fe (S J ) = 0if^j, (4) 

k 

where the Afe(Bj)'s are the eigenvalues of Bi with respect to the fixed ordering. 

The association scheme of particular interest to us arises as follows. Given a 
connected graph G = (V, E) with diameter d, we define \V\ x \V\ matrices A^ 
(k = 1, . . . , d) as follows: 

(fe) _ f 1 ifdist(i,j) = fc . , 

A ^ ~ \ else, {h] eV) > 

where dist(i,j) is the length of the shortest path from i to j. 

Note that is simply the adjacency matrix of G. Moreover, one clearly 

has 

i + Y j a^ = j. 

k=\ 

It is well-known that, for G = C n , the matrices A^ (k — 1, . . . ,d = [n/2\) 
together with A^ := I form an association scheme, since C n is a distance 
regular graph. 

It is shown in the Appendix to this paper, that for G — C n , the eigenvalues 
of the matrix A^ are: 

X m (A {k) ) = 2cos(27rmfc/n), m = 0, . . . , n - 1, k= 1, . . . , [(n - 1)/2J , 

and, if n is even, 

A„/ 2 (A (fc) )-cos(fc7r) = (-l) fc . 

In particular, we have 

A m (#')=A,(iW) fc,m=l,...,L(n-l)/2j. (5) 
Also note that 

\ m (A^) = \ n _ m (A^) 7 fc,m = l,...,L(n-l)/2j, (6) 
so that each matrix A^ (k = l,...,d) has only 1 + [n/2\ distinct eigenvalues. 
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Verifying the SDP relaxation ([3]) 

We now show that setting = A<- k > (k = l,...,d) gives a feasible solution 
of We only need to verify that 

j + ^cosf^W^O, i = l,...,d. 

k=l \ 71 / 

We will show this for odd n, the proof for even n being similar. 

Since the A^'s may be simultaneously diagonalized, the last LMI is the 
same as: 

d 

2 + YjA fc (A«)A J -(A«)>0, i,j = l,...,d, 

k=l 

and by using ([5]) this becomes: 

d 

2 + X k (A^)X k (A^) > 0, i.j 1 d. 



k=l 



Since Xq(A^) = 2 (i = l,...,d), and using (|4j) , one can easily verify that the 
last inequality holds. Indeed, one has 

d 

2 + ]TA fe (A«)A fc (#>>) 

fc=l 
1 11 — 1 

= 2 + -^A fc (AW)A fc (A«) (by©) 

fe=i 

1 1 11-1 

-X (A^)X (A^) + X k (A^)X k (A^) 



2 " ' v y 2 

2-2 + = if (i^j), by® 

2-2+iErd(A fe (^ (i) )) 2 >0 if = 

Thus we have established the following result. 

Theorem 1. The optimal value of the SDP problem ^ provides a lower bound 
on the optimal value TSP op t of the associated TSP instance. 



3 Relation of (JS) to an SDP relaxation of QAP 

An SDP relaxation of the QAP problem (JXJ) was introduced in [35] , and further 
studied for specially structured instances in [TJ. 

When applied to the QAP reformulation of TSP in ([2]), this SDP relaxation 
takes the form: 
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mm 
subject to 



\ trace(Ci <g> D)Y 



trace((J <g> (J - I))Y + ((J - J) <g> I)Y) = 
trace(F) — 2e T y = —n 
1 y T 
V Y 



>- o, y > o. 



(7) 



It is easy to verify that this is indeed a relaxation of problem @, by noting 
that setting Y = vec(X)vec(X) T and y = diag(Y) gives a feasible solution if 

ien„. 

In this section we will show that the optimal value of the SDP problem |7| 
actually equals the optimal value of the new SDP relaxation (|3|). The proof is 
via the technique of symmetry reduction. 



Symmetry reduction of the SDP problem ((7j) 

Consider the following form of a general semidefinite programming problem: 

p* := min { trace(^4o_X') : trace(^X) = bf., k = 1, ...,m}, (8) 

where the A.- L (i = 0, . . . , m) are given symmetric matrices. 

If we view ([7]) as an SDP problem in the form (JSj) , the data matrices of 
problem ([7]) are: 

( T W T W -e T \ ( 1 T \ 

\ \d®D )'\0 I®(J-I) + (J-I)®I, J'\ -e 21 j ' V )' 

(9) 

Definition 2. We define the automorphism group of a matrix Z S R kxk as 

aut(Z) = {P e U k : PZP T = Z}. 

Symmetry reduction of problem ([8|) is possible under the assumption that 
the multiplicative matrix group 

m 

G -=f] aut(A l ) 
i=a 

is non-trivial. We call Q the symmetry group of the SDP problem ([H]). 
For the matrices ©, the group Q is given by the matrices 

Q:={{ \ p ^) : (10) 

where T> n is the (permutation matrix representation of) the dihedral group of 
order n, i.e. the automorphism group of C n . 

The basic idea of symmetry reduction is given by the following result. 



G 



Theorem 2 (see e.g. [5]). If X is a feasible (resp. optimal) solution of the SDP 
problem ([8]) with symmetry group Q , then 



Pee 



is also a feasible (resp. optimal) solution of ([5]). 
Thus there exist optimal solutions in the set 

Ag := | 7^7 53 pT ^ p : x e K " X " | • 
I ' ' Pee J 

This set is called the centralizer ring (or commutant) of Q and it is a matrix 
^-algebra. For the group defined in (fT0|) , it is straightforward to verify that the 
centralizer ring is given by: 



Ac 



a x 
y C®Z 



G R, C = C T circulant. Z e 



x,y g 



(11) 

where x T = [xie T . . . x n e T ] and y T = [yie T . . . y n e T ] for some scalars Xi and 
(z = 1, . . . , n), where e G R™ is the all-ones vector, as before. 

Thus we may restrict the feasible set of problem to feasible solutions of 
the form (fTTjl . 

If we divide y and Y in ||7J into blocks: 



and 



V = 



Y 



(y {1) ) ■■■(y {n) ) 

y(n) ... y(i») 



yr(nl) _ _ , y(nn) 

where yW G R n and = Y^ T G M nx ™, then feasible solutions of satisfy 



/ i (y (1) ) 

y(n) 
y yW y( nl ) 



(y (n) ) T \ 
y(ln) 



y(™«) / 



>- 0. 



(12) 



Feasible solutions have the following additional structure (see [28] and The- 
orem 3.1 in [7]): 

• y W (i = 1, . . . , n) is a diagonal matrix; 

• yw) (i ^ j) is a matrix with zero diagonal; 
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• trace(jyW)) = 1 = 1, . . . , n); 

• Er=i^ fo) =e(2/ W ) T (i = l,...,n); 

• diag(r) = y. 

Since diag(F) = y for feasible solutions, we have j/W = diag(y( M )) (i = 
1, . . . , n). Moreover, since we may also assume the structure we have that 

y {t) =y % e (i = l,...,n), 

for some scalar values yi. This implies that the diagonal elements of yw all 
equal Since the diagonal elements of Y~(") sum to 1, we have £/, = 1/n and 
diag(y(")) = (l/n)e. Thus the condition: 



1 y 1 



y o 



reduces to 



Y-—JhO 



by the Shur complement theorem. This is equivalent to 

(I (g) Q*)Y{I ® Q) - \{I®Q*)J{I®Q) h 



where Q is the discrete Fourier transform matrix defined in (|25[) in the Appendix. 
Using the properties of the Kronecker product and of Q we get 



/ Q*Y^Q •• 


■ Q*Y^Q \ 




(k ■ 








-J® 






\ Q*Y^Q ■■ 


■ Q*y( ,i ")g j 




\ o • 





Recall that = ±1 and that we may assume yw) (j ^ j) to be symmetric 
circulant, say 

yW> (ijfei). 
fc=i 

where Cfc (k — l,...,d) forms a basis of the symmetric circulant matrices 
with zero diagonals (see the Appendix for the precise definition). Note that 
the nonnegativity of yW) is equivalent to x£ > (fc = 1, ...,d). Since 
trace(jy^- J - ) ) = 1 one has 



E 

fe=i 



(ij) 



2n 



Since J2t=x Y^ = e (j/^Tj — ^ J, one also has 



d n 



" " i 

EE^ = ~ J - 



fc=l ( = 1 
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By the definition of the Cfc's, this implies that 



E 



if 1 < k < [(n- 1)/2J 
if k = n/2 (n even). 



(13) 



Moreover, 



Q*YMQ=J2*rD k , (i^j), 



fc=i 



where Dk is the diagonal matrix with the eigenvalues (f2"r?)) of Ck on its diagonal. 
Thus the LMI becomes 



■ ELt4 ln) M 



J 



- J( 



V o 











^0. (14) 



The left hand side of this LMI is a block matrix with each block being a diagonal 
matrix. Thus this matrix has a chordal sparsity structure (n disjoint cliques of 
size n). We may now use the following lemma to obtain the system of LMFs 
©• 

Lemma 1 (cf. [Mj). Assume a nt x nt matrix has the block structure 

(£>(") ... £)(!») 
: '•■ : 
£)("!) . . . £)(«") 

where D™' £ St are diagonal = 1, . . . , n). Then M y if and only if: 
/ D {11) ••• D (ln) \ 

^ 7.7. 7.7. 



V 



■■■ Dl 

Applying the lemma to the LMI (fT4|) . and setting 

fe = l,...,Ln/2j, 



^ z = 1,. 



(15) 



yields the system of LMI's in ([3]). 

Thus we have established the following result. 



Theorem 3. The optimal values of the semidefinite programs ([3]) and J7]) are 
equal. 
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4 Relation of ([3]) to an SDP relaxation of Cvetkovic 
et al. 

We will now show that the new SDP relaxation ([3]) dominates an SDP relaxation 
(fl6|) due to Cvetkovic et al. [5] . This latter relaxation is based on the fact that 
the spectrum of the Hamiltonian circuit C n is known. In particular, the smallest 
eigenvalue of its Laplacian is zero and corresponds to the all ones eigenvector, 
while the second smallest eigenvalue equals 2 — 2 cos (— ). 
The relaxation takes the form: 

TSP op t > min i trace(L>X) 



subject to 



Xe = 2e, 

diag(X) = 0, 

< X < J, 

2J-X+(2-2cos(^))(J-Z) h 0. 



(16) 



Note that the matrix variable X corresponds to the adjacency matrix of the 
minimal length Hamiltonian circuit. 



Theorem 4. The SDP relaxation ((3]) dominates the relaxation (fT6|) . 
Proof. Assume that given X^ e S n (k = l,...,d) satisfy ©. Then, 

diagpf (1) ) = 0, 

while (Unj) and (TJ5]) imply 

xWe = 2e (k = 1, . . . , [(n - l)/2j), 

and X^l^e = e if n is even. In particular, one has X^e = 2e. It remains to 
show that 

21 - X<U + ^2 - 2 cos ^ ( J - I) h 0, 

which is the same as showing that 

21 - X™ + (2 ~ 2 cos (^) ) J2 X [k) >= 0, (17) 



k=l 



since 



fc=i 

We will show that the LMI p7[) may be obtained as a nonnegative aggregation 
of the LMI's 

d 

i + ^Tx {k) y 
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and 



I + ^cosf— W*) ^0 (i = l,...,d). 

k=l \ 71 / 

The matrix of coefficients of these LMFs is a (d + 1) x (d + 1) matrix, say A, 
with entries: 

' 2nij * 



A i:j = cos\— -J = Q,...,d). 

Since we may rewrite (| 1 T[) as 

2/ + f 1 - 2 cos f ) XV + f 2 - 2 cos f ) £ ^ 0, 



n I 



k=2 



we need to show that the linear system Ax = b has a nonnegative solution, 
where 



b := 



One may verify that, for n odd, the system Ax = b has a (unique) solution given 

by 

4 r d(l-cos(^)) if * = o 

„\ cos(^)-cos(^i) fari = l,...,d. 



Note that x is nonnegative, as it should be. If n is even, the solution is 

( t n— - 
4 ~ 



(ilil (i _ cos (1ZL)) if * = 



J., ■■; 



"I cos (f)^cos(^) fori = l,...,d-l 



□ 

In the section with numerical examples, we will present instances where the 
new SDP relaxation ((3J is strictly better than (|16p . 

5 Relation to the Held-Karp bound 

One of the best-known linear programming (LP) relaxations of TSP is the LP 
with sub-tour elimination constraints: 

TSP opt > min - trace(L>X) 

subject to 

Xe = 2, 
diag(X) = 0, . 

< X<J, ' (J,S) 

Z ie x,m X v ^ 2 VMZc{l,...,n}. 
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This LP relaxation dates back to 1954 and is due to Dantzig, Fulkerson and 
Johnson [6] . Its optimal value coincides with the LP bound of Held and Karp [15] 
(see e.g. Theorem 21.34 in [H]), and the optimal value of the LP is commonly 
known as the Held-Karp bound. 

The last constraints are called sub-tour elimination inequalities and model 
the fact that C n is 2-connected. Although there are exponentially many sub- 
tour elimination inequalities, it is well-known that the LP <| 1 8|) may be solved 
in polynomial time using the ellipsoid method; see e.g. Schrijver [53], §58.5. 

It was shown by Goemans and Rendl [T2] that this LP relaxation dominates 
the SDP relaxation (fT!)]) by Cvetkovic et al. [5]. The next theorem shows that 
the LP relaxation (|18[) does not dominate the new SDP relaxation ([3]), or vice 
versa. 



Theorem 5. The LP sub-tour elimination relaxation (|T8|) does not dominate 
the new SDP relaxation ([3]), or vice versa. 

Proof. Define the 8x8 symmetric matrix X as the weighted adjacency matrix 
of the graph shown in Figure [T] 

Edge weights: 
l 
2 




Figure 1: The weighted graph used in the proof of Theorem [5] 



The matrix X satisfies the sub-tour elimination inequalities, since the mini- 
mum cut in the graph in Figure [1] has weight 2. 

On the other hand, there does not exist a feasible solution of ([3]) that satisfies 
X^ = X, as may be shown using SDP duality theory. 

Conversely, in Section [7] we will provide examples where the optimal value 
of (|T8]) is strictly greater than the optimal value of ([3]) (see e.g. the instances 
grl7, gr24 and bays24 there). □ 

6 An LMI cut via the number of spanning trees 

In addition to the sub-tour elimination inequalities, there are several families 
of linear inequalities known for the TSP polytope; for a review, see Naddef [T5] 
and Schrijver [24 , Chapter 58. 
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Of particular interest to us is a valid nonlinear inequality that models the 
fact that C n has n distinct spanning trees. To introduce the inequality we require 
a general form of the matrix tree theorem; see e.g. Theorem VI. 29 in [Mj for a 
proof. 

Theorem 6 (Matrix tree theorem). Let a simple graph G — (V,E) be given 
and associate with each edge e <E E a real variable x e . Define the (generalized) 
Laplacian of G with respect to x as the \V\ x \V\ matrix with entries 

-x e ii{i,j} = e, 

else. 

Now all principal minors of L(G)(x) of order \V\ — 1 equal: 

T e(ET 

where the sum is over all distinct spanning trees T of G. 

In particular, if L(G)(x) is the usual Laplacian of a given graph, then x e = 1 
for all edges e of the graph, and expression ([19)) evaluates to the number of 
spanning trees in the graph. 

Thus if X corresponds to the approximation of the adjacency matrix of a 
minimum tour, then one may require that: 

det(27-X) 2:n2:n >n, (20) 

where X2-. n .2:n denotes the principle submatrix of X obtained by deleting the 
first row and column. 

The inequality (|20|) may be added to the above SDP relaxations (fT6|) and 
© (with X = XW), since the set 

{Z h : detZ > n} 

is LMI representable; see e.g. Nemirovski [19] . §3.2. 

We know from numerical examples that (|20[) is not implied by the relaxation 
of Cvetkovic et al. (I16p . but do not know any examples where it is violated by 
a feasible X^ of the new relaxation Nevertheless, we have been unable to 
show that ([HD (with X = X<U) is implied by ©. 

7 Numerical examples 

In Tabic [T] we give the lower bounds on some small TSPLir£] instances for 
the two SDP relaxations ((3]) and (fT6|) . as well as the LP relaxation with all 
sub-tour elimination constraints (fT5)) (the Held-Karp bound). These instances 

x http: //www. iwr.uni-heidelberg.de/groups/comopt/software/TSPLIB95/ 
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(19) 



have integer data, and the optimal values of the relaxations were rounded up to 
obtain the bounds in the table. 

The SDP problems were solved by the interior point software CSDP [2] using 
the Yalmip interface [T7] and Matlab 6.5, running on a PC with two 2.1 GHz 
dual-core processors and 2GB of memory. 



Problem 


SDP bound g]o$ 


SDP bound © (time) 


LP bound (JT8J) 


TSPopt 


grl7 


1810 


2007 (39s) 


2085 


2085 


gr21 


2707 


2707 (139s) 


2707 


2707 


gr24 


1230 


1271 (1046s) 


1272 


1272 


bays29 


1948 


2000 (2863s) 


2014 


2020 



Table 1: Lower bounds on some small TSPLIB instances from various convex 
relaxations. 



Note that the relaxation ([3]) can indeed be strictly better than {16]), as is 
clear from the grl7, bays24 and bays29 instances. Also, since the LP relaxation 
(fl~8|) gives better bounds than ([3]) for all four instances, it is worth recalling that 
this will not happen in general, by Theorem [5j 

The LMI cut from (|2H|) was already satisfied by the optimal solutions of p6[) 
and © for the four instances. 

A second set of test problems was generated by considering all facet defining 
inequalities for the TSP polytope on 8 nodes; see [31 for a description of these 
inequalities, as well as the SMAPO project web sitqj. 

The facet defining inequalities are of the form ^ trace(-DX) > RHS where 
J}£cS„ has nonnegative integer entries and RHS is an integer. From each in- 
equality, we form a symmetric TSP instance with distance matrix D. Thus the 
optimal value of the TSP instance is the value RHS. In Table[2]we give the opti- 
mal values of the LP relaxation (fT8| (i.e. the Held-Karp bound), the SDP relax- 
ation of Cvetkovic et al. (flB|) , and the new SDP relaxation © for these instances, 
as well as the right-hand-side RHS of each inequality h tra,ce(DX) > RHS. For 
n = 8, there are 24 classes of facet-defining inequalities. The members of each 
class are equal modulo a permutation of the nodes, and we need therefore only 
consider one representative per class. The first three classes of inequalities are 
sub-tour elimination inequalities. 

The numbering of the instances in Table [21 coincides with the numbering of 
the classes of facet defining inequalities on the SMAPO project web site. 

The new SDP bound © is only stronger than the Held-Karp bound (fl"8j) 
for the instances 16, 21 and 23 in Table [2j and for the instances 1, 5, 9, 12, 14, 
17, 18, 19 and 22 the two bounds coincide. For the remaining 18 instances the 
Held-Karp bound is better than the SDP bound ©. However, if the bounds are 
rounded up, the SDP bound ([3]) is still better for the instances 16, 21 and 23, 
whereas the two (rounded) bounds are equal for all the other instances. Adding 

2 http : //www . iwr . uni-heidelberg . de/ groups/comopt/sof tware/SMAPO/tsp/ 



14 



I n cc]u a,l i t y 
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Table 2: Results for instances onn = 8 cities, constructed from the facet denning 
inequalities. 



the LMI cut from (f2T)]) did not change the optimal values of the SDP relaxations 
HU) or §5§ for any of the instances. 

For n — 9, there are 192 classes of facet defining inequalities of the TSP 
polytope [4]. Here the SDP bound ((3|) is better than the Held-Karp bound 
for 23 out of the 192 associated TSP instances. Similar to the n = 8 case, 
when rounding up, the rounded SDP bound remains better in all 23 cases and 
coincides with the rounded Held-Karp bound in all the remaining cases. 

8 Concluding remarks 

Wolsey [57] showed that the optimal value of the LP relaxation (|18[) is at least 
2/3 the length of an optimal tour for metric TSP (see also [5S]). An interesting 
question is whether a similar result may be proved for the new SDP relaxation 

Finally, the computational perspectives of the SDP relaxation ([3]) are some- 
what limited due to its size. However, since it provides a new polynomial-time 
convex approximation of TSP with a rich mathematical structure, it is our hope 
that it may lead to a renewed interest in improving approximation results for 
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metric TSP. 
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Appendix: Circulant matrices 

Our discussion of circulant matrices is condensed from the review paper by Gray 
A circulant matrix has the form 



C 



C Ci c 2 • • • 
c n -i c ci 

Cn-l Co Ci 



Cl 



Cn-l 



Cl 

Cn-l C 



Thus the entries satisfy the relation 

Cu = 



^3 \3~i) mod n' 



The matrix C has eigenvalues 



\ m (C) = c + J2 c k e- 27rV ~ lmk/n , m = 0, . . . , n - 1. 



k=l 



(21) 



(22) 



If C is symmetric with n odd, this reduces to 

(n-l)/2 

A m (C) = Co+ 2cfe cos(27rm/c/n), m = 0, . . . , n — 1, (23) 

fc=i 

and when n is even we have 

n/2-l 

A m (C)=c + 2cfe cos(27rmfc/?i) + c„/ 2 cos(m7r), m = 0, . . . , n — 1. (24) 
fc=i 

The circulant matrices form a commutative matrix *-algebra, as do the sym- 
metric circulant matrices. In particular, all circulant matrices share a set of 
eigenvectors, given by the columns of the discrete Fourier transform matrix: 



Qi 



.— J_ e -2n^Tij/r 



i,j = 0, ...,n- 1. 



(25) 



1G 



One has Q*Q = I, and Q*CQ is a diagonal matrix for any circulant matrix C. 
Also note that Q*e — y/ne. 

We may define a basis C^-°\ . . . , CL"/ 2 J for the symmetric circulant matrices 
as follows: to obtain we set c, = Cn-i = 1 in (f2"Tj) and all other Cj's to zero. 

(We set Co = 21 and also multiply C„/ 2 by 2 if n is even). 

By (|2"51) and (f2"4")l , the eigenvalues of these basis matrices are: 

A m (C (fc) ) = 2cos(27rmfc/n), m = 0, . . . , n - 1, k = 0, . . . , ln/2\ . (26) 
Also note that 

A m (C^) = A„_ m (C«), m = 1, . . . , [n/2\, k = 0, . . . , [n/2\ 
so that each matrix has only 1 + [n/2\ distinct eigenvalues. 
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